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ABSTRACT 

In a protoplanetary disk, a combination of thermal and non-thermal desorption processes regulate 
where volatiles are liberated from icy grain mantles into the gas phase. Non-thermal desorption should 
result in volatile-enriched gas in disk-regions where complete freeze-out is otherwise expected. We 
present ALMA observations of the disk around the young star IM Lup in 1.4 mm continuum, C^®0 2- 
1, H^^CO+ 3-2 and DCO“*' 3-2 emission at ~0"5 resolution. The images of these dust and gas tracers 
are clearly resolved. The DCO"'" line exhibits a striking pair of concentric rings of emission that peak 
at radii of ^0"6 and 2" (^90 and 300 AU, respectively). Based on disk chemistry model comparison, 
the inner DCO+ ring is associated with the balance of CO freeze-out and thermal desorption due to a 
radial decrease in disk temperature. The outer DCO’*' ring is explained by non-thermal desorption of 
CO ice in the low-column-density outer disk, repopulating the disk midplane with cold CO gas. The 
CO gas then reacts with abundant H 2 D+ to form the observed DCO^ outer ring. These observations 
demonstrate that spatially resolved DCO+ emission can be used to trace otherwise hidden cold gas 
reservoirs in the outmost disk regions, opening a new window onto their chemistry and kinematics. 
Subject headings: astrochemistry; line: profiles; molecular processes; techniques: imaging spec¬ 
troscopy; protoplanetary disks; circumstellar matter; ISM: molecules; radio lines: 
ISM 


1. INTRODUCTION 

Planets form through the assembly of refractory dust, 
volatile ice, and gas in disks around young stars. The 
spatial distribution of volatile abundances helps deter¬ 
mine the compositions of nascent planets that form at 
different disk radii. The separation of these volatiles into 
their gas and ice phases is set by the balance of adsorp¬ 
tion (“freeze-out” of gas onto solids), and desorption back 
to the gas phase. Desorption results from both thermal 
and non-thermal processes that can be active in different 


regions of protoplanetary disks. 

The thermal desorption (sublimation) rate for a 
volatile is set by its energy barrier for ice desorption in 
concert with the grain temperature. Common volatiles 
like H 2 O, CO 2 , and CO are characterized by very dif¬ 
ferent desorption barriers. For a typical disk with mid¬ 
plane temperatures that increase closer to the host star, 
this range of desorption energy barriers will introduce a 
series of distinct locations - “snow lines” - where spe¬ 
cific volatiles are liberated from ices into the gas phase. 
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The highly volatile species CO and N 2 have the low¬ 
est thermal desorption barriers, and so have snow lines 
at large disk radii where the disk is cold, ~30-100AU 
(|Oberg et al.l l2011all . Irradiation heating maintains a 
thermal inversion in the disk atmosphere, where temper- 
atures increase wit h height (z) above the midplane (e.g. 
iCalvet et alJll99Jl . These higher surface temperatures 
mean that snow lines are really two dimensional snow 
surfaces whose boundaries extend to larger r for larger 
z. 

In addition to thermal desorption, there are analogous 
non-thermal processes driven by interactions between ice 
and high-energy photons and particl es, or the release of 
chemical energy (jCarrod et al.l 120071 1 . Ice photodesorp¬ 
tion by ultraviolet (UV) radiation has garnered special 
interest because of the strong UV fields produced by 
young stars (e.g., iBergin et al.l l2003ll and the high UV 
desorption yields measured in laboratory experiments 
(iWestlev et al.l 119951: lOberg et al.1 l2007t iFavolle et al.1 
l2011l : lChen et al.ll2014ll . In the inner disk, UV photodes¬ 
orption can result in lower scale heights for snow surfaces 
compared to the on e s expected from thermal desorption 
(IWoitke et al.l 120091 : iHogerheiide et al.l 120111 : lOka et aP 
l2012fl . At the low densities of the outer disk, UV radi¬ 
ation may penetrate to the cold midplane to desorb CO 
ice well o utside the expected locat ion of the thermal CO 
snowline (|Willacv fc Langeill2000ll . 

DCO'*' is expected to form efficiently in cold regions 
where CO gas is still present in the gas-phase, whether 
thermally or non-thermally desorbed, through deuteron 
transfer from the low-temperature (T<50 K) product 
H 2 D+ to CO. Therefore, the expectation is that DCO"^ 
emission peaks close to the CO snow surface, i.e. at 
the coldest disk locations where CO is still abundant in 
the gas-phase if only thermal desorption is considered. 
Several studies have found an increasing 000"*" abun¬ 
dance with radius a scribe d to this process (lOi et al.ll2nn8l 
iTeague et al.ll2015^ . and iMathews et ^ (1201311 resolved 
a ring of DCO^ in the disk around the Herbig Ae star HD 
163296, which seems closely related to the CO snow line. 
Their proposed connection between the DCO'*" emission 
morphology and the CO snowline is challenged by new 
disk chemistry models, where significant DCO'*' can be 
produced at temperatures a.bove 30 K throug h pathways 
that do not involve H 2 D“'' (|Favre et al.ll201^ . However, 
this complication does not eliminate the fact that DCO+ 
requires gas-phase CO to form, nor does it remove the 
expectation that DCO+ is formed in excess if H 2 D+ is 
abundant (i.e., if the gas is cold). 

In this study, we present ALMA Cycle 2 observations 
of the IM Lup disk, located at a distance of 155 pc 
(|Lombardi et al.ll200^ . in 1.4 mm continuum, DCO+ 3- 
2, C^®0 2-1 and H^^CO'*" 3-2 line emission (§2). §3 
presents the observed emission morphology of DCO'*' and 
the other species. In §4 we use a disk chemistry model, 
applied to a generic T Tauri disk structure, to interpret 
the observations. We highlight the role of non-thermally 
desorbed CO to explain the presence of DCO+ emission 
in the outer regions of the disk in the discussion in §5 
and summarize our findings in §6. 


2. OBSERVATIONS AND DATA REDUCTION 


Observations of the IM Lup disk were ac¬ 
quired with the Atacama Large Millime¬ 
ter/Submillimeter Array (ALMA; project code 
ADS/JAO.ALMA#2013.1.00226.S) on 2014 July 6 
with 31 antennas (targeting DCO’*' and C^®0) and on 
2014 July 17 with 32 antennas (targeting H^^CO’*'). 
The antenna separations spanned baselines of 15-650 m. 
The total on source integration time was ~21 minutes 
in each execution. The nearby quasars J1534—3526 and 
J1427—4206 were used for gain and bandpass calibra¬ 
tion, respectively. Titan was observed to calibrate the 
absolute amplitude scale. 

On the July 6 execution, the correlator was configured 
to process 13 distinct spectral windows (SPWs). The 
spectral resolution Jz/ and bandwidth for 12 of the SPW 
were 61 kHz and 59 MHz, respectively; the remaining 
SPW had a coarser resolution with a bandwidth of 469 
MHz. The DCO"*" 3-2 line was located in SPWl, and 
C^®0 2-1 in SPW6 (Table [T|). For the July 7 execution, 
the correlator was configured to process 14 SPWs in a 
different part of the Band 6 spectrum. All had the same 
spectral resolution, 61 kHz; twelve SPWs had a band¬ 
width of 59 MHz, and the remaining two had a band¬ 
width of 117 MHz. The H^^CO'*' 3-2 line (Table [IJ was 
located in SPWl. 

The visibility data were calibrated by ALMA staff. 
The calibration plots were inspected, deemed adequate, 
and no changes to the delivered calibration were made. 
Each individual SPW was further phase and amplitude 
self-calibrated in CASA 4.2.2. In each SPW, the con¬ 
tinuum was subtracted using line-free channels. The 
fully calibrated visibilities were Fourier inverted and 
CLEANed using the Briggs weighting scheme: the ro¬ 
bust parameter was set 0.5 for the continuum and 1.0 for 
the spectral lines. A CLEAN mask was manually gener¬ 
ated based on the bright C^®0 2-1 line, and then applied 
to the H^^CO'*' and DCO+ lines. The CLEANed maps 
were restored with a synthesized beam with FWHM di¬ 
mensions of ^ 0'.'7 X 0'.'5. The 1.4 mm continuum has an 
RMS noise level of ^0.24 mJy beam“^ (dynamic range 
limited). The integrated flux density is 181 ±18 mJy 
(which includes a 10% flux calibration uncertainty), con¬ 
sistent with previous SMA measurements (jOberg et al.l 
I2011bll . The resulting RMS noise level in 0.2 km s ^ 
binned channels is 4-6 mJy. 

3. OBSERVATIONAL RESULTS 

Eigure[T] shows the continuum and integrated line emis¬ 
sion maps of C^®0 2-1, H^^CO’*' 3-2, and DCO"*" 3-2 
from the IM Lup disk. The continuum is centrally peaked 
and spatially resolved, extending out to a radius of ~2" 
(300 AU), comparable to previous lower resolution obser¬ 
vations (IPinte et al.ll2008HPanic et al.ll2009HOberg et al.l 
I2011bfl . The peak flux density at the disk center is 66±7 
mJy beam“^ at 1.4 mm (C^®0 and DCO+ spectral set¬ 
ting) and 78±8 mJy beam“^ at 1.1 mm (H^^CO+ spec¬ 
tral setting), corresponding to a brightness temperature 
of ~6 K. 

Integrated line emission maps (Fig. [TJ top panel) were 
constructed by summing up emission from individually 
CLEANed 0.2 km s“^-wide channels, employing a 2cr 
threshold to maximize the SNR. They are shown together 
with the continuum 2cr contour. The middle panel of 
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TABLE 1 

Molecular line data 


Molecular Line 

Rest freq. 

GHz 

Logic (Aij) 

Eu 

K 

beam (PA) 

" X " (°) 

Integrated flux 
mJy km/s 

DCO+ J=3-2 

216.1126 

-2.62 

20.7 

0'.'65 X 0l'48 (-71°) 

450 ± 45 

Ci®0 J=2-l 

219.5604 

-6.22 

15.8 

0'.'68 X 0l'47 (-60°) 

1460 ± 150 

H13cO+ J=3-2 

260.2553 

-2.87 

25.0 

0'.'68 X 0"59 (89°) 

484 ± 48 



3 2 1 0 -1 -2 -3 

Aa ["] 
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Fig. 1.— Continuum and integrated emission maps of 2—1, 3-2 and DCO"*" 3-2 toward IM Lup (top row), integrated 

emission in two velocity bins (middle row) and deprojected azimuthally averaged profiles (bottom row), top row: The contours in the 
continuum map are 2cr+[2,4,8,16,...]cr with the continuum rms a = 0.24 mjy. The 2cr continuum contour is also shown on top of the 
integrated line emission maps. The integrated line emission in units of mJy km/s beam“^ is shown in gray scale (see individual color bars 
in each panel). The continuum peak is marked by a red cross and the synthesized beam is displayed in the bottom left of each panel. 
middle row: The integrated emission in two velocity bins around the source velocity has been scaled by 0.3 for and 0.7 for HCO'^ 

to enable all three line maps to be shown on the same scale (color bar to the right of the DCO'^ panel), bottom row: The azimuthally 
averaged molecular emission profiles have been labeled to highlight the lack of molecular emission at the continuum maximum, and the 
clear DCO*^ double -ring structure. 


Fig [T] presents integrated line emission maps of the same 
lines split up into two velocity bins around the source 
velocity to visualize the disk rotation. The integrated 
emission in these panels have been scaled to a constant 
emission maximum to show the increasing importance of 
extended emission for DCO+ compared to H^^CO+, and 
for compared to C^®0 . To visualize the radial 

structure of these tracers more clearly and to mitigate 
any bias introduced by the applied clipping threshold, 
we also produced moment-0 maps without any clipping, 
deprojected and azimuthally averaged them into radial 
brightness profiles (Fig. [U bottom panels). We assumed 
an inclinatio n of 49° and a position angle of -35° f or these 
calculations (jPinte et al.ll2nn8HPanic et al.ll2n09[l . 


The integrated emission maps and average radial pro¬ 
files reveal two significant morphological features. First, 
all three molecular species show a clear central depres¬ 
sion. The C^®0 and DCO^ emission each peak at a 
radius of 0'.'6, or ^90 AU, and the line peaks at 

0"8 (~125AU). Second, DCO+ exhibits a large second 
ring at 2", near the outer edge of the continuum emis¬ 
sion. No such ring is observed in C^®0, or the 

continuum. However, the H^^CO+ emission “plateaus” 
around the same radius and C^®0 exhibits a subtle slope 
change. 

4. DISK CHEMISTRY MODELING 
4.1. Model Description 
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To interpret the observed C^®0, H^^CO+, and DCO+ 
radial structures, we explored the rele vant predictions of 
the di sk chemistry mo d el pre sented bv iFuruva &: Aikawal 
(j2Q14f l: lAikawa et all (HoH), updated with deuterium 
chemistry that considers nuclear spin (ortho/para) chem¬ 
istry of light species. This model adopts a generic 
steady, axisymmetric Keplerian disk around a T Tauri 
star with stellar mass, radius, and effective temperature 
of M* = 0.5 Mq, ii* = 2 Rq, and = 4000 K, resp^ - 
tively, and disk mass 1.7x 10“^ Mq (jNomura et al.ll2?)0^ . 
The stellar UV and X-ray luminosities are set to 10^^ erg 
s“^ and 10^° erg s“^, respectively, b ased on the observed 
spectrum of the you ng star TW Hya (|Herczeg et akllW)^ 
iKastner et al.lf2?i?)^ . The disk model assumes that the 
grains have grown up to 1 mm throughout the disk, re¬ 
sulting in a factor of 10 decrease in grain surface area 
compared to ISM dust. The gas temperature, dust tem¬ 
perature and density distributions of the disk are cal¬ 
culated self-consisten tly, considering vario us heating and 
cooling mechanisms (|Nomnra et al.ll2007fl . 

Compared to IM Lup, this model employs a fainter 
star, which results in a colder disk, and also a lower disk 
mass. The differences in temperature and density struc¬ 
ture clearly precludes any quantitative comparison be¬ 
tween the model and data. However, even though the 
model has not been tuned to the parameters of IM Lup, 
qualitative comparisons still provide valuable insight into 
the origins of the observed molecular emission structure, 
as long as the model results are robust to the model pa¬ 
rameter choices. To evaluate the latter we performed a 
small parameter study where we varied the CO binding 
energy, the disk density and temperature structure, and 
the cosmic ionization rate. As described below the main 
results are not sensitive to these parameters. We there¬ 
fore conclude that the predictions from the generic disk 
model are useful for interpreting the observed emission 
pattern of the specific disk around IM Lup. 

The chemical composition of the model is calcu¬ 
lated by integrating a time-dependent system of gas 
and grain surface rate equation s for 300 kyr, assumin g 
the physical structure is static (iHaseeawa et al.l [199211 . 
Th e gas-grain reacti o n net work is based on the work 
by iCarrod fc Herbstl (|2n06[l . supplemented with calcu¬ 
lations th at account for hig h-temperature gas phase 
reacti ons (iHarad a et al.l 1201(1(1 . X-ray mediated chem¬ 
istry (iFuruva et al.lI2013L and re f erenc es therein), deu¬ 
terium chemistry (lAikawa et al.l I2012L and references 
therein), and nuclear s pin state chemistry of H^, Hj~ 
and t heir isotopologues (IHueo et 1111120091: iCoutens et al.l 
[Ml . Species containing chlorine, phosphorus, or more 
than four carbon atoms were excluded for computa- 
tional expediency. Elem ental abundances are taken from 
lAikawa fc Herbstl (120011 ). The initial abundances are ob¬ 
tained by calculating t he molecular evoluti on of a star¬ 
forming core, following lAikawa et al.l (I2012D . 

Importantly for this study, the binding energy of CO 
in the generic disk model is set to 1150 K, its mea - 
sured value for H 2 O ice surfaces (|Collings et al.l (Mil , 
which results in a CO midplane freeze-out temperature 
of ~26 K. If CO is instead freezing out on a pure CO 
ice, the binding energy will be lower. We explored 
the effects of such a scenario by running a model with 
the CO binding energy set to 1000 K. Non-thermal 


desorption i s included through stoch astic heating by 
cosmic-rays ([Hasegawa fc Herbs^ll993D , photodesorption 
(lOberg et al.1 l2009al lbl: iFavolle e t al.l 120 111: Che n et al.l 
120141 ). and reactive desorption ( Garrod et al.l 2007f|. The 
CO p hotodesorption yield is set to 7 x 10“^ (jFavolle et al.l 
1201,11 ). UV photodesorption and cosmic ray heating are 
the most important non-thermal desorption mechanisms. 
Their relative contributions depend on a combination of 
adopted yields, cosmic ray ionization rate, dust proper¬ 
ties, and UV field, which varies across the disk. 

The main target of this study, DCO^, forms through 
three different pathways in the model 

1. H2D+-hCO^ H2-kDCO+, 

2. HC0+-kD-4DC0+-kH 

3. CH 2 D+-I-O—^DC0+-|-H2 and other pathways re¬ 
lated to CH 2 D+. 

The D and CH 2 D+ pathways regulate DCO'*' formation 
at elevated temperatures, i.e. in the inner (R<10 AU) 
disk and in the disk atmosphere at all radii. The D path¬ 
way is generall y the more importa nt of the two in the 
present model. iFavre et al.l (|2015D . on the other hand, 
argued that the CH 2 D'^ path is more important, proba¬ 
bly due to a different disk ionization structure compared 
to our model. Exterior to 10 AU, warm DCO'*' is a mi¬ 
nor contribution to the total DCO'*' column density. In¬ 
stead, the H 2 D+ pathway dominates the DCO'*' produc¬ 
tion. H 2 D+ forms efficiently at low temperatures due to 
the small (<230 K) zero-point energy difference between 
Hg and H 2 D+ (the exact energy difference depends on 
the spin states of H 2 and H 2 D+). 

4.2. Model Results and Comparison to Observations 

Figure[2]shows the model disk temperature and density 
structure. One particular difference worth highlighting 
between this generic model and IM Lup is in the disk 
midplane temperature profile. In the model, 20 K is 
reached at a radius of ~40 AU, while a previous analysis 
of the IM Lup disk suggests that 20 K is reached at a 
larger radius (^40- 100 AU), depende nt in detail on the 
model assumptions (|Pinte et al.ll2008fP1 . As a result, CO 
freeze-out in this model probably starts at smaller radii 
compared to the IM Lup disk. 

These model calculations do not include C and O 
isotope-specific chemistry. To make comparisons with 
the observations of isotopes, we have scaled the calcu¬ 
lated CO and HCO+ densities by factors of 500 and 70, 
respectively, adopting the standard ISM ^®0/^®0 and 
12 c/13c, estimate the C^®0 and H^^CO'*' densities. 
We expect that this is a reasonable approximation in the 
disk midplane (in the disk atmosphere, isotopologue spe- 
cihc photodissociation rat es will change the re lative CO 
isotopologue abundances (jMiotello et al.ll2?)T^ ). 

Figure [2] (middle row) shows the radial and vertical 
gas-phase abundance profiles for the species that are the 
focus of this study, i.e. CO, HCO"*" and DCO'*', with the 
first two abundances scaled to correspond to the observed 
isotopologues. All three molecules exhibit an abundance 
depression in the disk midplane exterior to 25 AU due to 

^ We have corrected for the assumed distance of 190 pc 
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Fig. 2.— Density, gas and dust grain temperatures in our generic disk model (top row) and the CO/500 (~C^®0), HCO’*"/70 (~H^®CO^) 
and DCO+ abundances with respect to H nuclei after SOOkyr (middle row). The density contours mark nH=10^ and 10® cm“® (upper 
left panel). The temperature contours mark 20, 30, 50 and 100 K (upper middle and right panels). The bottom row shows the calculated 
radial emission profiles of the observed lines C^®0 2-1, H^®CO+ 3-2 and DCO+ 3-2. 


freeze-out of CO onto grain surfaces (i.e., the midplane 
CO snow line). Between 25 and 50 AU, the CO abun¬ 
dance decreases by seven orders of magnitude, the HCO+ 
and DCO'*" abundance by five orders of magnitude. 

Exterior to 90 AU, the CO, HCO’*' and DCO"*" gas- 
phase midplane abundances begin to increase. With 
decreasing density, the freeze-out rate decreases, which 
changes the CO desorption/freeze-out balance and there¬ 
fore increases the gas/solid ratio. The accompanying de¬ 
crease in column density results in increased UV penetra¬ 
tion and therefore CO ice UV photodesorption, further 
increasing the CO gas/solid ratio. The result is a sub¬ 
stantial amount of CO gas in the cold, outer disk mid¬ 
plane, corresponding to a C^®0 abundance of 10“^^ 
with respect to nn- This is six orders of magnitude above 
the abundance minimum found at 50 AU, but still small 
compared to the total C^®0 abundance of ~ 10“^, i.e. 
most CO (isotopologues) are present in the form of ice in 
the outer disk midplane. By contrast, a C^®0 abundance 


of ^ 10“^ is reached in the warm molecular disk layer, 
ZjR ^ 0.3 — 0.5, and this is the main reservoir of C^®0 
gas in the outer disk. 

1300+ and DCO'*' are also present in the warm molec¬ 
ular disk layer, but this layer is no longer the dominant 
reservoir. In the case of DCO'*', the midplane and molec¬ 
ular layer abundances are comparable. The relatively 
high DCO"*" abundance in the midplane, at ZjR < 0.2, 
is due to the increased efhciency of DCO’*' production 
at low temperatures. Because of an increasing gas den¬ 
sity toward the disk midplane, the vertically integrated 
DCO'*' column density is dominated by the midplane 
DCO'*' in the outer disk. HCO'*' presents intermediate 
trends between CO and DCO'*'. 

We calculated disk emission spectral data cubes us- 
ing the 3D Monte Carlo rad iative transfer code LIME 
(|Brinch fc Hogerheiidell2010l ) from the model disk den¬ 
sity, temperature and molecular abundance structures, 
for the observed C^®0, H^®CO’*' and DCO'*' lines. Figure 
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[H (bottom row) shows the radial emission profiles (arbi¬ 
trary units) of the C^®0, and DCO'*' obtained 

from the model spectral-image data cubes, integrated in 
velocity space and azimuthally averaged like the data. 
The C^®0 emission profile is centrally peaked. It shows 
a radial gradient small change at ^25 AU at the onset 
of CO freeze-out. The effect of CO gas in the midplane 
beyond 90 AU on the CO emission profile is subtle, and 
in practice unobservable, since the overall column den¬ 
sity and thus emission at these radii is dominated by 
the abundant gas-phase CO in the warmer atmosphere 
layers. 

The HCO+ and DCO+ emission profiles decrease 
rapidly at radii of 25-40 AU, due to the freeze-out of 
gas-phase CO (their formation reservoir). The HCO+ 
emission then presents a “plateau” out to 100 AU, at¬ 
tributed to Cosmic Ray and UV photodesorption of CO 
back into the gas phase. The effect of the midplane CO 
non-thermal desorption on the DCO+ emission is, as ex¬ 
pected, more dramatic. The return of CO to the gas at 
large radii, where low temperatures drives ann efficient 
deuterium fractionation chemistry, results in a second 
maximum (outer ring) peaking at 200 AU. DCO'*' also 
presents a small peak toward the central star, caused by 
warm deuterium fractionation chemistry, and an inner 
ring close to the CO midplane snow line. Unlike for CO 
emission, the presence of cold CO in the outer midplane 
has a considerable effect on HCO+ and DCO'*' emission, 
which should be observable. 

Based on a small grid of models the shapes of the pre¬ 
dicted CO, HCO’*' and DCO+ radial profiles are not sen¬ 
sitive to specific model assumptions. In particular the 
two concentric DCO^ rings are reproduced when chang¬ 
ing the CO binding energy from 1150 to 1000 K, i.e. from 
a water to a CO dominated ice environment, and when 
the CR ionization rate is reduced to 10“^° s“^. In the lat¬ 
ter case UV photodesorption still maintains cold CO gas 
in the outer disk, albeit at a lower level compared to when 
a high cosmic ray ionization rate is assumed. That is, the 
outer ring is still present when the cosmic ray ionization 
is attenuated, but it may be less pronounced. The molec¬ 
ular emission trends in Fig thus appears to be robust 
as long as some non-thermal desorption process is active 
in the disk midplane. However, the absolute DCO’*' col¬ 
umn density and predicted emission profile in the outer 
disk do depend on model parameters, and whether the 
predicted DCO+ double-ring will be observable may thus 
vary from disk to disk. 

Qualitatively, the observed and predicted CO, HCO+ 
and DCO+ radial profiles agree well. In particular, the 
model reproduces several key features of the IM Lup ob¬ 
servations: the centrally peaked C^®0 emission profile, 
the H^^CO+ shoulder in the outer disk, and the DCO+ 
inner ring and outer ring. However, the model shows no 
signs of a central depression for C^®0, and only a very 
small depression for H^^CO+. For DCO+, the model 
double-ring is accompanied by an emission peak at the 
disk center, which is also not observed. 

5. DISCUSSION 

5.1. The DCO'^ Outer Disk Ring 

The most surprising result from the ALMA observa¬ 
tions is the clear detection of a concentric pair of emis¬ 


sion rings in the DCO'*' 3-2 spectral line, peaking at radii 
of ~90 and 300 AU. Based on the midplane temperature 
profile for the IM L up disk that was i nferred from prev i- 
ous measurements (iPinte et al.ll2008HPanic et al.ll2009ll . 
the inner ring may coincide with the onset of CO freeze- 
out. The relationship between the CO midplane snowline 
and the inner DCO'*" ring location is, however, compli¬ 
cated by multiple DCO^ formation pathways in the inner 
disk. 

The outer DCO'*' ring requires the return of some CO 
into the gas-phase at large disk radii. This can be ac¬ 
complished by either non-thermal ice desorption or a ra¬ 
dial thermal inversion in the low density outer region. 
Of the two scenarios we favor the former, since it nat¬ 
urally explains the large DCO'*' enhancement and the 
much smaller H^^CO'*" enhancement at large disk radii. 
This could be observationally checked by observing ad¬ 
ditional DCO"*" and H^^CO'*' lines that constrain the gas 
temperature. 

The location of the outer DCO'*' ring coincides with 
the disappearance of mm-grain emission (Fig. [1]). It may 
also coincide with a steepening of the gas-density profile 
(jPanic et al.ll200^ . The radial profile of smaller grains, 
which carry most of the opacity at UV wavelengths, is 
not known independently, but it seems reasonable to as¬ 
sume that it too is steepening at the same radius. The 
DCO+ outer ring then likely coincides with a disk lo¬ 
cation where the UV radiation penetration depth is in¬ 
creasing rapidly with radius. The connection between 
the DCO'*' outer ring and UV radiation is supported by 
chemical modeling. An outer DCO'*' ring is produced 
when CO is non-thermally desorbed from icy grains in 
the cold, low-density outer disk. In the model UV pho- 
todesorption and CR heating are both important, but the 
coincidence between the ring and the loss of dust opacity 
suggests that in the case of IM Lup, UV photodesorp¬ 
tion is the most important desorption mechanism in the 
outer disk. Once in the gas phase, CO can be chemically 
converted into DCO"^ following deuteron transfers from 
abundant reservoirs of H 2 D+ and other deuterated H;^ 
isotopologues (which are enhanced at such low temper¬ 
atures). This effect is expected to be generally present 
in disks and especially pronounced in disks with signif¬ 
icant grain growth, which reduces UV opacity. In our 
model we assume a uniform grain growth throughout the 
disk, parameterized as a reduction in total grain surface 
area by a factor of 10. Further grain growth is likely 
present in older disks, and it may therefore be illustra¬ 
tive to compare the DCO'*' distribution in disks of dif¬ 
ferent ages and with different grain properties as esti¬ 
mated from continu um observations (|Perez et al.ll201^ 
[Andrews et al.l[2f)^ . 

Based on the discussion above and the seemingly ro¬ 
bust model results in M.21 it is tempting to conclude that 
concentric DCO'*" rings should be commonly observed in 
protoplanetary disks. It is important to realize however, 
that the expected in tensity of the secon d DCO'*' ring 
is model dependent. iTeague et al.1 (j2015l l. for example, 
finds a DCO'*' enhancement in the outer disk, but it is 
small compared to the f indings in our model. The same 
is true in the model bv iWillac^ (I2007D . In general, the 
DCO'*' abundance in the outer disk will depend on the 
intensity and penetration depths of radiation fields (UV, 
cosmic rays and X-rays), and the DCO'*' ring column 
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density will depend on the total amount of cold molec¬ 
ular gas exposed to this radiation. The intensity of the 
DCO“*' ring will also depend on the excitation conditions 
that characterize the outer disk midplane, especially its 
density and temperature. It may very well be that 000"^ 
rings, while common, are only observable in a small sub¬ 
set of disks where disk density, temperature and radia¬ 
tion structures conspire to produce sufficient amounts of 
DCO+ flux at large disk radii. Still the detection of one 
such outer ring in the IM Lup disk demonstrates that 
non-thermal desorption in general, and UV photodes¬ 
orption of CO in particular, is important in setting the 
chemical structure of protoplanetary disks. 

5.2. Origins of the Central Depression in Several 
Tracers 

A lingering puzzle in the observations presented here 
is the central depression in the intensities of all the ob¬ 
served spectral lines (but notably not in the continuum). 
We consider four potential explanations for these depres¬ 
sions: (1) a reduction in the gas surface densities; (2) 
line-of-sight absorption by foreground cloud material; (3) 
unexpected chemistry; and (4) high continuum optical 
depths that prevent some of the line emission from es¬ 
caping the disk. 

The first possibility seems unlikely, since there is no 
analogous signature in the continuum, and we are un¬ 
aware of a physical mechanism that would systematically 
enhance the mm dust-to-gas ratio in the inner disk. The 
second scenario, molecular cloud contamination, should 
produce deeper depressions (absorption) for the more op¬ 
tically thick spectral lines. Inspection of CO data toward 
the same source show that CO exhibits less visible de¬ 
pressions than C^®0, ruling out this scenario. 

Chemistry can produce central depressions in HCO'*" 
and DCO+ following desorption of H 2 O and other 
molecules with high proton affinity in the inner disk, 
resulting in efficient proton transfer from HCO^ and 
DCO’*' to H 2 O. Chemistry may thus contribute to the 
H^^CO"*" and DCO+ holes. The C^®0 central depression 
is not predicted by our or any other publish e d disk chem¬ 
istry model (e.g. iWillacy fc WoodsI 120091 iWalsh et al.l 
I2OIOI : iSemenov fc Wiebel l2011[ l. though the conversion 
efficiency of CO into other species across the disk is an 
active topic of research fe.g. lFavre et al.l[20I3[l . Based on 
the results presented so far, it seems unlikely, however, 
that chemistry alone can explain the entire observed ef¬ 
fect. 

The last proposed explanation is that the continuum 
optical depths in the disk center are high enough to 
block some of the line emission from the midplane. The 
deeper depressions observed for more optically thin trac¬ 
ers, which originate closer to the midplane, and for lines 
emitted at higher frequencies (especially H^^CO+) sup¬ 
ports this scenario. But this scenario has one significant 
problem: the peak continuum intensity at the disk cen¬ 
ter corresponds to a brightness temperature of only 6 K, 
much lower than expected for optically thick emission 
at these radii. The only way of reconciling this expla¬ 
nation with the data is through significant beam dilu¬ 


tion. If the optically thick continuum emission is con¬ 
centrated in small scale features (occupying an area ^5 x 
smaller than the synthesized beam), then this scenario 
may work. That may seem unlikely, but the recent obser- 
vations of the HL Tau disk show ing small scale structure 
(jALMA Partnership et al1l2ni5h indicate that it remains 
a possibility. This hypothesis could be tested directly by 
very high spatial resolution continuum measurements, or 
less directly by observations of CO and HCO+ isotopo- 
logue lines at 3 mm, where the dust opacity should be 
reduced. 

6. SUMMARY AND CONCLUSIONS 

ALMA observations of the massive disk around IM Lup 
have revealed a spectacular DCO+ double-ringed struc¬ 
ture. The inner ring appears connected with the ther¬ 
mal CO snowline, but dust opacity and multiple DCO^ 
formation pathways complicate the interpretation. The 
second, outer DCO’*' emission ring at 300 AU is read¬ 
ily reproduced by non-thermal CO ice desorption in the 
low-density, outer disk. The ring coincides with a rapid 
decrease in dust opacity, which is likely accompanied by 
increased UV penetration toward the disk midplane. We 
hypothesize that the DCO+ outer ring is caused by effi¬ 
cient UV photodesorption of CO ice at this radius. 

Without the presence of non-thermal desorption this 
outer disk region would be completely devoid of CO gas 
and thus of CO-gas-mediated chemistry. Even at its most 
efficient, non-thermal desorption in disk mid planes is not 
expected to maintain a large fraction of CO in the gas- 
phase, < 0.1% of the total CO reservoir. The presented 
observations and models show, however, that the release 
of such trace amounts of CO into the gas-phase can still 
produce an observable signature through its effects on 
DCO'*" production. This DCO+ can be used as a probe 
of the kinematics and ionization in the outer disk mid¬ 
plane, opening up a new window on the structures of 
protoplanetary disks. 
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